Doping dependence of spin fluctuations and electron correlations in iron pnictides 
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Doping dependence of the spin fluctuations and the electron correlations in the eiTective five-band 
Hubbard model for iron pnictides is investigated using the fluctuation-exchange approximation. For 
a moderate hole doping, we flnd a dominant low-energy spin excitation at Q = (tt, 0), which becomes 
critical at low temperature. The low-energy spin excitations in the heavily hole-doped region are 
characterized by weak Q dependence. The electron doping leads to an appearance of a pseudogap in 
spin-excitation spectrum. Correspondingly, the NMR-l/Ti relaxation rate is strongly enhanced on 
the hole-doped side and suppressed on the electron-doped side of the phase diagram. This behavior 
can be to large extent understood by systematic changes of the Fermi-surface topology. 

PACS numbers: 74.20.-z; 74.25. Jb; 74.70.Xa; 76.60.-k; 78.70.Nx 



I. INTRODUCTION 



The recent discovery of iron-pnictide superconductors 
with high transition temperatures (Tc)r^ some over 
50K)^ has provoked an intense research. The appearance 
of high- Tc superconductivity in a close proximity to the 
anti-fcrromagnctic (AF) phase, induced by a carrier dop- 
ing and an applied pressure, reminds of the cuprates and 
the heavy-fermion superconductors. The sign-reversing 
s^-wave state with the pairing glue of AF fluctuations 
has been suggested as an explanation of the observed 
superconductivity.-"— The NMR— "— and the neutron- 
scattering experiments^^"— revealed in detail the evolu- 
tion of magnetic excitations with the carrier doping, in 
particular the presence of the strong AF fluctuations in 
the hole-doped systems and opening of a (pseudo)gap 
for the electron doping. The uniform susceptibility has 
shown puzzling linear increase over a broad temperature 
range, which does not match with the two common lim- 
its, T- independent for weakly interacting electrons and 
1/T for local moments. The effect of electron correlations 
has been observed and quantified by means of the anglc- 
rcsolved photoemission spectroscopy (ARPES)r^ the op- 
tical spectroscopy?^ and de Haas-van Alphen (dHvA) 
experimentsj^i^ The accumulated knowledge suggests 
that understanding the changes of spin dynamics and 
electronic correlations with the carrier doping provides 
the key to the superconducting-paring mechanism. 

Previously, one of us (H.L) studied the doping and 
temperature dependence of the single-particle spectra 
and the NMR-l/Ti, applying the fluctuation-exchange 
(FLEX) approximation to an effective five-band Hubbard 
model, obtained for the band structure of LaFeAsOii 
This early study succeeded in predicting the enhance- 
ment of the AF spin fiuctuations on the hole-doped side, 
but also revealed a problem with double counting the in- 
teraction effects in the multiband systems. Straightfor- 



ward addition of the FLEX sclf-encnrgy leads to redis- 
tribution of orbital occupancies from their LDA values 
accompanied by drastic changes in the Fermi surface and 
the spin fluctuations in contrast to experimental observa- 
tions. In the preceding papersj^i^ we have investigated 
this point in detail and proposed a way around, subtrac- 
tion of a static part of the self-energy, together with phys- 
ical argument supporting this ad hoc procedure. With 
this modiflcation, we now can apply the FLEX approxi- 
mation over wide ranges of carrier doping and tempera- 
ture. 

Here we employ the modified FLEX scheme^ to study 
the effect of varying electron concentration on spin- 
fiuctuations in the five-band model of iron pnictides 
and compare it with the experimentally observed trends 
across the pnictide series. In order to single out the ef- 
fect of carrier doping we do not construct precise tight- 
binding models for each individual compound, but use 
the parameters obtained for LaFeAsO<^ The calculations 
capture the common features in the series of iron-pnictide 
superconductors derived from LaFcAsO and BaFe2As2. 



II. MODEL AND COMPUTATIONAL DETAILS 

Wc start with a five-band Hubbard model in the un- 
folded two-dimensional Brillouin zone (BZ). The kinetic 
term 
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comes from the tight-binding model in LaFeAsO, and 
the hopping integrals appear in TABLE II of Ref. HI. 
Only the on-site Coulomb interaction is considered with 
the common parametrization: the intraorbital Coulomb 
U, the interorbital Coulomb U' , the Hund coupling J, 
and the pair- hopping J'. In the FLEX approximation. 
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the normal Green's functions Qgrn{k,iuJn) for orbitals £ 
and m are self-consistently determined from the fohowing 
equations, 

Gimik) = g'^^ik) + J2 Gu'{k)^e'm'{k)g,n'm{k), (2a) 

t'm' 

^hn{k) = X! X! ^U',mm'iq)gtm'{k - q) , (2b) 



(a) 



(b) 



q Cm' 



'{q) = [uu~2Un~Unx''{q)Un 

+ ^U^nqP^ + lu'r{q)U' 



(2c) 



Here the bare vertices U^''^ = Ufi =F U^^ with {U^i)u,u = 

U , (C/ti)«,™™ — {Uti)hn,im ~ J, iU^l)im,mi = J', 

iU-f^)u,mm = U' - J, {U^-l)em,im = J ~ U' fOT £ ^ ITl, 

and the susceptibihties in the spin sector and the charge 
sector are given by 



r{q)=x''{q)+x'{q)U'r{q), 
nq)^x'{q)~x"{q)U'r{q), 



(3a) 
(3b) 



and 



[X°(l)]u'.,mm' = - E ^^"(^ + l)Grn'cik). (4) 



With the straightforward appHcation of FLEX, we en- 
counter the double counting problem mentioned in the 
Introduction. Therefore, as in Ref. IH, we subtract the 
Lu = part of the self-energy in order to eliminate 
the unwanted redistribution of the charge between or- 
bitals. In Ref. m we argued that this mimics the effect 
of the Hartree part of the electron-electron interaction, 
which does not appear in our low-energy effective model. 
The correction term Y,f^{k,u! = 0) was calculated at 
T = 23K (Ref. and used unchanged also at higher 
temperatures so that the low-temperature Fermi surface 
matches the LDA one, but at higher temperature small 
modifications of the Fermi surface are allowedi^ Once 
the correction term is fixed, this procedure equals the 
standard FLEX approximation for a Hamiltonian with 
the kinetic part, = - T,f^{k, 0) instead of 

Through this paper, we take 64 x 64 meshes in the 
unfolded BZ and 1024 Matsubara frequencies, and U = 
1.20eV and J = 0.25eV as the interactions with U ~ 
U' + 2J and J' = J . The retarded quantities gf^{k,uj), 
S^(fc,a;) and xf(QiW) are obtained by the numerical 
analytic continuation with use of the Pade approxima- 
tion. 



III. RESULTS AND DISCUSSION 

A. Single-particle spectra and Fermi surface 

The Fermi Surfaces (FS) formed by the renormalizcd 
bands for fillings^ corresponding to heavy, n = 5.52, 
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FIG. 1: (Color online) Contour plots of the Green's function, 
-5^Jme«(fe,i7rr)/7r, at T = 4meV~ 46K for n = 5.52 (a), 
5.92 (b), and 6.16 (c) in the unfolded BZ. Deep red corre- 
sponds to the Fermi surface. The Fermi surface around M 
point is the electron sheet at n = 5.92 and 6.16, but the hole 
sheet at n = 5.52. The weight around F' point becomes small 
at lower temperatures, since the band edge locates below the 
Fermi level. 



and moderate, n = 5.92, hole-doping and heavy electron- 
doping, n = 6.16 are shown in Fig. [TJ The broadening 
comes from the imaginary part of the self-energy, arising 
mainly from the magnetic scattering. The n = 5.92 FS 
consists of hole sheets around the F and the F' points and 
an electron sheet centered at the M point nested with the 
well known vectors Q = (tt, 0) or (0, tt). Electron doping 
leads, as expected, to shrinking of the F {F') hole sheets 
and expansion of the M point electron surface, as shown 
in Fig. [TJ:. At n = 6.16 the F' surface is reduced to 
almost a point, which leads to a T dependence discussed 
below. The n = 5.52 FS, similar to that of KFe2As2, has 
large hole sheets around the F {F') points and a small 
hole pocket around the M point. 

Figs. [2^ andj^b show the orbital- resolved spectral den- 
sities pi{uj) = — X]fclni^«(fci i^)/''' for = 5.52 and 
6.16. Like for the non-interacting bands, the Fermi level 
E-p = falls to a vicinity of a spectral peak, more so 
for the hole doped systems. The states on the Fermi 
surface arise predominantly from d,^z/yz and d,^2_y2 or- 
bitals. Their relative weight varies strongly with doping 
and while on the hole-doped side the d^z/yz contribution 
dominates over d^ja-ya, the reverse holds for the electron- 
doping. 

Figs. [2j: and [2ji show the T-dependent changes of the 
spectral functions in the vicinity of the Fermi level. A re- 
markable suppression of the dx-2_yi density at Ey upon 
cooling can be seen for n = 6.16, when the Fermi level 
is located close to a step-like van Hove singularity at the 
top of the F' band. The decreasing temperature leads 
not only to reduced quasi-particle damping, but also to 
a downward band shift and the corresponding shrinking 
of the F' pocket. Such a variation of the spectral density 
has been observed in the photoemission spectra;^ and 
numerical study of Ref. 0. The shrinking of FS has been 
also observed in LaFePO (Ref.lli) and BaFe2(Asi_a;Pa;)2 
(Ref. 113) ■ General arguments for correlation-induced 
band shifts due to coupling of an asymmetric electronic 
band to a bosonic mode, e.g. spin fluctuation, were given 
by Ortenzi et al^ Applying their reasoning together with 
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around n ~ 5.92 where nesting, and thus strongly k- 
dcpendent effects, dominate the physics. In Fig. we 
show the doping dependence of the mass enhancement 



factors for each orbital, 



ImZ'ff (fc, inT) 



(5) 



in the vicinity of the F and M sheets of the Fermi surface. 
Wc observe an overall trend of increasing mass enhance- 
ment from 2 — 3 on the electron-doped side to over 4 on 
the hole-doped side, which can be attributed to a growing 
spectral density at Ep (see Fig. [5^). In agreement with 
this trend, a remarkably large mass enhancement fac- 
tor is observed in dHvA on KFe2As2»^ As shown below, 
this can be because low-energy spin fluctuations become 
featureless, Q-independent, similar to the heavy fermion 
systems. The cusps in for dj.2_y2 orbital at n 5.92 
result from strong stripe-type AF spin fluctuations due 
to the FS nesting. 



FIG. 2: (Color online) The partial density of states pt{i^) 
(states/eV) at n = 5.52 (a) and n = 6.16 (b) at 46K. Temper- 
ature dependence for d^^/y^ and d^2_y2 orbitals at n = 5.52 
(c) and n — 6.16 (d). Shown is remarkable suppression of 
Pt{0) for d^2_y2 orbital at n = 6.16. 
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FIG. 3: (Color online) (a) The self-energy for d^2_y2 orbital 
at the r' point at 46K in units of eV. (b) Mass enhancement 
factor zl. for each orbital at fe = (267r/32,0) around the M 
point (solid) and fc = (67r/32, 0) around the F point (dotted). 



the increase of spin- fluctuation density at low-T (Fig. U) , 
indeed, leads to the observed shift of the F' band. 

The main effect of electronic correlations in itinerant 
systems is an enhancement of the quasiparticle mass. 
In Fig. [3^, we use typical self-energies to demonstrate 
that the scattering processes contained in FLEX lead 
to quasiparticle renormalization over a rather broad en- 
ergy range. In fact, we find a fair agreement between 
the FLEX self-energy and the self-energy obtained with 
LDA-fDMFT for the five band modcl4S This applies also 
to the mass enhancement factors, except for the region 



B. Spin-excitation spectra 

Next, we discuss the structure of the spin excitation, 
Imxf(Q,a;), at n ~ 5.52, 5.92, and 6.16, displayed 
in Fig. 01 In the left panel we show contour plots of 
Imxf'{Q,uj) along the high symmetry directions. The 
right panel shows the variation of Imxf (Q, w) at fixed Q 
with temperature. 

For a moderate hole doping, n = 5.92, the strong low- 
energy fluctuations (stronger than in the undoped case), 
arising from FS nesting, are located around Q = (tt, 0) 
(see Fig. HJd). They exhibit strong enhancement at low T 
indicating an incipient AF instability (Fig. HJ^). Further 
hole doping destroys the nesting between F {F') and M 
sheets of FS, and increases the density of states at the 
Fermi level as shown in Fig. [5^. This leads to rather fea- 
tureless, Q-independent, and only weakly T-dependent 
structure of the low-energy spin fluctuations as shown in 
Figs. 2^, HJi, and|3^. A broad hump around Q — (7r,0), 
which develops at low F, is due to scattering between 
the F hole sheets and the M pocket, composed of d^z/yz 
orbitals, in Fig. [1^. 

In contrast, the electron doping results in suppression 
of the low-energy spin fluctuations. For n = 6.16 the 
(tt, 0) spin fluctuations are remarkably suppressed, and 
a gap opens in the spin excitation spectrum at low T as 
shown in Figs.^j; andUf. Recently, such a gap- like behav- 
ior with heavy electron doping has been observed by the 
inelastic neutron-scattering experiment This behav- 
ior is linked to the corresponding changes of the Fermi- 
surfacc topology in Fig. [TJ:. As the electron doping leads 
to the shrinking of the F {F') hole sheets and expansion 
of the M electron sheet, the particle-hole excitation at 
Q = (tt, 0) requires a finite energy, and the correspond- 
ing spectral weight moves to higher energies. The lead- 
ing scattering channel becomes dominated by scattering 
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FIG. 4: (Color online) Imaginary part of xfiQ:'^) (A'b/sV) along the high-symmetry line at n = 5.52 (a), n = 5.92 (b), and 
n = 6.16 (c). Temperature dependence at Q = (7i",0) and (7r,7r/2) at n = 5.52 (d, g), n = 5.92 (e, h), and n — 6.16 (f, i). 
The system at n = 5.92 locates at around the stripe-type AF critical point. At n — 5.52, shown is Q-independent low-energy 
excitation, while the gap- like behavior at n = 6.16. 



between different electron (M) sheets, and the domi- 
nant low-energy spin fluctuation moves to Q = (tt, 7r/2), 
as shown in Fig. |3J. Previous calculations showed that 
the leading superconducting instability changes from s±- 
wave to dx2_y2-wa,ve around this doping.— This points 
to the correlation between the structure of the spin fluc- 
tuations and the superconducting pairing symmetry. 

C. NMR l/Ti and uniform susceptibility 

Next, we discuss the local spin response mea- 
sured in terms of the NMR-l/Ti relaxation rate. 
Fig. [S^ shows the temperature dependence of l/TiT = 
^qlmxf (q,a;)/w|(^^o for various dopings with the 
hyperfine-coupling constant set to unity. The overall 
trend of 1 /TiT growing with the hole doping at high tem- 
peratures follows from increasing pe(0). At lower temper- 
atures, l/TiT is suppressed on the electron-doped side, 
reflecting the reduction of p«(0), as shown in Fig. [2ji, 
and the corresponding opening of a spin gap (Fig.|3t and 
llf). Such suppression of 1/TiT with electron doping, 
which is consistent with the previous works,^^^ has been 
observed in LaFeAs(0i_a;F2;) (Ref. fla - flSr). LaFeAsOi_y 
(Ref. Hi) and Ba(Fei_^Co^)2As2 (RefllJ. A remark- 
able growth of l/TiT in n = 5.92 and 5.80 cases comes 



from the enhanced AF spin fluctuation shown in Figs.HjD 
and 21:. Although such enhancement is suppressed with 
further the hole doping, a slight increase can be ob- 
served also for n — 5.52. originating from the broad 
hump structure around Q ~ (tt, 0) in Figs. HJi and 
UJi. The trends in l/TiT are consistent with a remark- 
able low-T increase in Ban.7 9Kn .28Fe2As2 (Ref. [isl ) and 
Bao.6Ko.4Fe2As2 (Ref. dalM, and a relatively small 
increase in KFe2As2 (Ref. [2l|[23 ) at low temperatures. 

Finally, we present T-dependence of the uniform sus- 
ceptibility xf(0,0) in Fig. [5}d. Roughly speaking, it ex- 
hibits an increase with temperature independent of the 
carrier doping except for n = 6.00 case. On the hole- 
doped side, the tendency is accompanied by the devel- 
opment of the AF fluctuation, in qualitative agreement 
with a scenario studied in Ref. |43|. On the electron- 
doped side, it is consistent with the suppression of the 
NMR-l/TiT and pe{Q), which is not accompanied by a 
remarkable development of the magnetic fluctuation. Ex- 
perimental observations shown clear a temperature de- 
pendence on the electron-doped side in LaFeAsOi-^-Fj; 
and Ba(Fei_j:Co2:)2As2,— i^*^ and a constant behav- 
ior on the hole-doped side in Bai_a:Ka;Fe2As2f^'^ Thus, 
we need further investigations for complete understand- 
ing of the uniform susceptibility. 
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FIG. 5: (Color online) (a) NMR l/TiT as a function of tem- 
perature. The vertical line is a logarithmic scale. Shown are 
remarkable enhancement on the hole-doped side and suppres- 
sion on the electron-doped side, (b) The uniform susceptibil- 
ity x?(0, 0) (/Xe/eV), which overall decreases upon cooling.— 



IV. CONCLUSIONS 

In conclusion, we have investigated the normal-state 
properties in the effective five-band Hubbard model 
for the iron-pnictides using the FLEX approximation. 
We have obtained a variety of trends in the spin dy- 
namics and the electron correlations with the carrier 
doping, which qualitatively agree with the overall fea- 
tures observed in the (Ba,K)Fe2As2, Ba(Fe,Co)2As2 and 
LaFcAsO systems, in particular, the gap-like feature in 
(tt, 0) spin-excitation spectrum in the heavily electron- 
doped case, corresponding to the Ba(Fei_2;Co2;)2As2 
with large x, and the weak Q dependence of the spin 
fluctuation and the large mass enhancement in the heav- 
ily hole-doped case, corresponding to the end material 
KFe2As2. We find that the changes in the Fermi-surface 
topology are the main driving force behind the observed 
trends. 
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